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ABSTRACT 


More than 17,000 measurements were taken on 
over 1000 specimens of two species of brachio- 
pods, Mediospirifer audaculus and Athyris spirifer- 
oides from the Middle Devonian Hamilton Group 
of New York State. Statistical analyses were per- 
formed in order to test patterns and processes of 
stasis in the morphology of species over several 
million years. Measurements were taken both to 
quantify shell shape and organismal ontogeny. 
Specimens were partitioned by their occurrence in 
one of many environments and stratigraphic ho- 
rizons. Stasis could not be refuted, as neither spe- 
cies showed substantial morphological departures 
between the oldest and youngest samples of the 
Hamilton Group (a roughly 5-million-year inter- 
val). However, oscillations in morphology, allo- 
metric heterochrony, and rates of change of single 
features were discovered in both taxa. 


In addition, one of the mechanisms typically 
held to be responsible for morphological stasis, 
stabilizing selection in conjunction with habitat 
persistence, does not appear to be the sole cause 
of stasis in these taxa. For the two species we stud- 
ied, groups of organisms occurring in a single en- 
vironment undergo substantial morphological 
change through time, suggesting that a stable en- 
vironment is not producing morphological stabil- 
ity. However, the net sum of changes through time 
in all the environments these species occur in is 
essentially zero. Therefore, stasis appears to be 
partly a property of the organization of species. 
Species are organized into different environmental 
populations or demes. Different “environmental 
populations” may evolve, but they will typically 
do so in several different “directions” in mor- 
phospace, generally producing no net change. As 
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long as a species occurs in several different envi- 
ronments it will be resistant to change. The dif- 
ference between the morphology of species in dif- 
ferent environments over the whole interval of the 
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Hamilton Group is also nil, thereby ruling out any 
major role that ecophenotypic effects could play 
in the patterns recognized herein. 


INTRODUCTION 


The pattern and rate of morphological 
change in evolutionary species lineages over 
geological time has been a hotly debated topic 
in paleontology and biology over the last two 
decades. The traditional, or syntheticist, 
model of evolutionary change within species, 
propounded in the work of Haldane (1931), 
Simpson (1953), Mayr (1963), and others, 
and referred to as phyletic gradualism, differs 
significantly from the work of the punctuated 
equilibrists, most importantly Eldredge and 
Gould (1972), Gould and Eldredge (1977), 
Stanley (1979), etc., in the relative impor- 
tance that the two theories ascribe to the 
prevalence of morphological stasis in species 
over long periods of time. 

Publication of the punctuated equilibria 
hypothesis touched off many endeavors to 
test the rate of change in species through time. 
Some studies have indicated that over time 
stasis is the rule (Williamson, 1981; Stanley 
and Yang, 1987; Gould, 1988; Barnosky, 
1990; Lich, 1990; etc.). Others have suggest- 
ed that gradual trends do occur (Gingerich, 
1976, 1985; Sheldon, 1987). 

This analysis takes as its starting point the 
search for a gradualistic pattern of irrevers- 
ible and orthogenetic morphological change 
in order to test the hypothesis of stasis. It also 
assesses the role that the environment plays 
in mediating stasis and change through time. 
To attain these goals, two of the most abun- 
dant species in the Middle Devonian Ham- 
ilton Group of New York State, the brachio- 
pods Mediospirifer audaculus (Conrad) and 
Athyris spiriferoides (Eaton) were studied. 
These species are found in almost all of the 
stratigraphic units of the Hamilton Group, 
are plentiful throughout the section, and oc- 
cur in a range of paleoenvironments. 

The Hamilton fauna, and the strata of the 
Hamilton Group, have been the subject of 
many detailed stratigraphic and paleoenvi- 
ronmental studies (e.g., Brett et al., 1990: 
Grasso, 1986; Miller, 1986; Savarese et al., 
1986), and these studies allow fossil speci- 


mens collected from the Hamilton Group to 
be assigned a relatively precise age and to be 
placed in a distinct paleoenvironment or bio- 
facies. The intensely studied geological sys- 
tem of the Hamilton Group presents a good 
opportunity to look at the relationship be- 
tween morphological change, stratigraphic 
history, and environment in individual spe- 
cies lineages. 

Punctuated equilibria and the competing 
phyletic gradualist model each make clear 
predictions about the anticipated patterns of 
morphological change in a single species 
through time. The former predicts that spe- 
cies will not change morphologically for most 
of their existence, while the latter argues that 
morphology will gradually, and in a mostly 
unidirectional fashion, change over long pe- 
riods of time (Eldredge and Gould, 1972; 
Gould and Eldredge, 1977). In both theories 
the abiotic environment is ascribed an im- 
portant role as a mechanism either driving 
or restricting evolutionary change in a species 
over time. 

The environment has been related to pat- 
terns we see in species morphology in at least 
three different ways: 1) it maintains stasis via 
stabilizing selection, 2) it produces orthoge- 
netic change by initiating a constant selective 
regime (assumes organisms are not yet adapt- 
ed to their environment), and 3) environment 
and morphology are decoupled such that 
morphology undergoes a random walk re- 
gardless of environment (unless environment 
is severe enough to induce death). The first 
explanation has typically been touted as one 
of the main causes of stasis (e.g., Lande, 1980; 
Charlesworth et al., 1982). (Alberch, 1982; 
Maynard Smith et al., 1985; and Williamson, 
1987 provided excellent discussions of an- 
other chiefly cited cause of stasis, develop- 
mental constraints/species integrity.) 

If stabilizing selection played a prominent 
role in canalizing the morphology of species, 
we would predict that species in constant en- 
vironments would not change, if similar en- 
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vironments exert similar selective regimes 
over time. When a brachiopod species from 
the Hamilton Group is found in the same 
biofacies or paleoenvironment at two differ- 
ent times, it is held to be occurring in the 
same environment. If stabilizing selection is 
a valid explanation of stasis then as long as 
the members of a species can persist in their 
preferred environment, they will not change. 
Thus, in similar environments, populations 
should not change over time. 

We performed a series of multivariate sta- 
tistical analyses on morphometric data from 
two of the most abundant brachiopod species 
restricted in age to the Hamilton Group in 
order to evaluate the prevalence of stasis. We 
also conducted analyses on individual mor- 
phological variables, to quantify their rates 
of change through time. In addition, we de- 
termined the amount of morphological change 
in species lineages through time across their 
entire sampled distributions and within sin- 
gle environments, to assess the role that sta- 
bilizing selection within a single biofacies or 
paleoenvironment may play in mediating 
stasis. The information on phenotypic evo- 
lution within an individual paleoenviron- 
ment/biofacies was considered in light of 
techniques that relate changes in phenotype 
to minimum selective mortality regimes. 


MATERIALS AND METHODS 


GEOLOGICAL AND PALEOECOLOGICAL 
SETTING CONSIDERED 


The geological setting we used to test the 
hypothesis of stasis in species lineages is the 
Hamilton Group, a package of rocks from 
New York State that contains around 350 
species, and spans 5 million years of the Mid- 
dle Devonian. Some of these taxa appear to 
have Armorican and North African affinities 
(Bailey, 1983; Eldredge, 1985) and their ap- 
pearance in North America is hypothesized 
to have been driven by the post-Eifelian 
(Middle Devonian) reintensification of the 
Acadian Orogeny (Brett, 1986; Cooper et al., 
1942; Ettensohn, 1985). From its inception, 
the Hamilton Group and its fauna are asso- 
ciated with a drastic increase in the amount 
of sedimentation from mountains to the east 
(Brett, 1986; Cooper et al., 1942; Ettensohn, 
1985). These mountains formed during the 
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Acadian Orogeny, which was caused either 
by the collision between Armorica, Lauren- 
tia, and the intervening Traveler terrane 
(Kent, 1985); by oblique convergence of Ava- 
lonia and Laurentia (Ettensohn, 1985); or by 
collision between Armorica-Iberia and Ava- 
lonia with subsequent effects on Laurentia 
(Soper et al., 1992). 

The erosional products of this orogenic 
event imparted a pattern of increasing ter- 
restrialization through time and in a west-to- 
east gradient in the Hamilton Group rocks. 
The thickest sections are found in eastern New 
York at the Catskill front and measure 1000 
meters. Near the westernmost extremity of 
New York State, the Hamilton Group is 100 
m thick (Brett, 1986). The bulk of the samples 
surveyed are from central and western New 
York, as it is there that the marine facies of 
the Hamilton Group are prominently devel- 
oped. Further east, nonmarine environments 
are found. 

The section spans latest Eifelian to middle 
Givetian and deposition commenced at about 
380 Ma and ended around 375 Ma on the 
basis of Rb-Sr dates (Brett, 1986; Brett and 
Baird, 1986). Deposition took place in the 
northern arm of the large, tropical epeiric sea 
that covered much of eastern North America 
(Brett et al., 1986) (fig. 1). A stratigraphic 
section for the Hamilton Group is shown in 
figure 2, along with bars indicating the geo- 
logical units from which specimens of M. au- 
daculus and A. spiriferoides could be ob- 
tained. Samples were grouped according to 
their occurrence in well-circumscribed 
chrono- and lithostratigraphic units. 

Fossil assemblages from the Hamilton 
Group can be classified into approximately 
eight generalized biofacies or paleoenviron- 
ments which were defined by factor analytic 
techniques using the presence and abundance 
of associated species as well as other criteria 
given in Brett et al. (1990), Grasso (1986), 
Miller (1986), and Savarese et al. (1986). 
There is cyclic change in the type of fossil 
assemblage present, which approximately 
follows a pattern of lithofacies (rock type) 
cyclicity (Savarese et al., 1986), and this cy- 
clicity records the upslope and downslope 
migration of assemblages tracking habitat 
shifts and changes in water depth (Brett et 
al., 1990). This, along with other evidence, 
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Fig. 1. Schematic diagram of the depositional 
environment of the Middle Devonian Hamilton 
Group of New York State, with a generalized sketch 
of the paleogeography showing the paleolatitude 
of the North American craton (Laurentia) at this 
time (380 million years ago). Modified from Brett 
et al. (1986). 


indicates that these biofacies are sound pa- 
leoenvironmental indicators. (For more de- 
tails on the designation of biofacies see Brett 
et al. [1990].) Locality information was used 
to assign specimens to one of the Hamilton 
paleoenvironments/biofacies discussed in 
Brett et al. (1990) (see fig. 3). The biofacies 
considered were: the Athyris-Mediospirifer 
biofacies, abbreviated as (B) (where the 
Athyris-Mediospirifer biofacies is taken as 
equivalent to the Athyris biofacies of Brett et 
al. [1990]), the Athyris-Mediospirifer transi- 
tional to Pseudoatrypa biofacies (C), the 
Pseudoatrypa biofacies (D), the Tropidolep- 
tus biofacies (E), the Pentamerella-Helio- 
phyllum biofacies (F), and the Mucrospirifer 
transitional to Ambocoelia biofacies (G). The 
number of specimens of Mediospirifer au- 
daculus and Athyris spiriferoides measured in 
each paleoenvironment/biofacies is given in 
figure 3. (There are slightly fewer specimens 
considered in the analysis by communities 
because the exact biofacies could not always 
be determined.) 
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MEASUREMENTS TAKEN AND 
ADJUSTMENTS FOR SIZE 


Data were gathered from each brachiopod 
Shell with a digital caliper accurate to 0.01 
mm. All measurements were rounded off to 
the nearest 0.1 mm. Measurements were taken 
using both homologous landmarks and max- 
imum length/width measurements. 

For M. audaculus, six measurements were 
taken on the pedicle valve at its terminal on- 
togenetic stage. This valve was chosen over 
the brachial valve because it contains greater 
amounts of utilizable information and is more 
resistant to preservational deformation. The 
measurements taken are illustrated in fig- 
ure 4. 

401 specimens were measured, and these 
specimens were divided into samples repre- 
senting 10 different stratigraphic horizons that 
spanned basal to terminal Hamilton strata. 
To avoid any biases, all available specimens 
were measured. The position of the samples 
in the section is shown by the darkened bars 
in figure 2, with the number of specimens for 
each stratigraphic horizon indicated in the 
figure. 

Data were adjusted to remove differences 
between horizons governed by differences in 
size using two different methodologies in or- 
der to cross-verify their efficacy. This was 
done to minimize the potential impact size 
has on the analysis because of the impact 
ecophenotypic variation has on size, and to 
correct for possible taphonomic biases relat- 
ing to size. In the first method, data were log 
transformed. In the second, magnitudes of 
the 1 X 6 specimen vector were transformed 
using the variables H, JAW, and MAXF__SW 
which in regressions against each other had 
slopes equal to 1.0, indicating isometric bi- 
variate slopes (see Somers, 1989; Sundberg, 
1989). The correction factor used the square 
root of the sum of the squares of these vari- 
ables. After generating derived data in which 
all variables were adjusted for size using the 
correction factor, the derived data were plot- 
ted against the size factor. For all variables, 
no association was found between the cor- 
rection factor and the scores of variables, in- 
dicating that for all variables size had been 
largely adjusted for. Both methods produced 
congruent results, but for the purposes of 
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et al. (1942), showing the major lithostratigraphic units, the horizons treated for each species, and the 
number of specimens measured for each of these species. The number of specimens of Athyris spiriferoides 
from the following horizons are: G = 45, H = 52, I = 54, K = 81, L = 35, M = 44, N = 35. 
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Fig. 3. Diagrammatic depiction of the paleoenvironments/biofacies from the Hamilton Group that 
were sampled in this analysis (modified from Brett et al. [1990]). The top of the diagram represents the 
shallowest environments and the right side of the diagram represents increasingly turbid water conditions. 
The number of specimens of M. audaculus analyzed for each biofacies are: from the Athyris/Mediospirifer 
biofacies (abbreviated B in the text) 84, the Athyris/Mediospirifer transitional to Pseudoatrypa biofacies 
(abbreviated C in the text) 77, the Pseudoatrypa biofacies (D) 152, the Tropidoleptus and Tropidoleptus 
transitional to Pseudoatrypa biofacies (E) 35, the Pentamerella/Heliophyllum and Pentamerella/Helio- 
phyllum transitional to Pseudoatrypa biofacies (F) 33. The number of specimens of A. spiriferoides 
measured are: (B) 162, (C) 45, (D) 243, (E) 76, (F) 40, and the Mucrospirifer transitional to Ambocoelia 


biofacies (G) 35. 


brevity and clarity only the results using the 
second method are presented here. 

Four measurements were made on 614 
specimens of Athyris spiriferoides (see fig. 5) 
from 18 different stratigraphic horizons. All 
measurements were taken on the pedicle valve 
because it is less subject to preservational de- 
formation. Data were scaled to reduce the 
divergence between horizons caused by size 
differences by employing the two method- 
ologies used for M. audaculus. In the vector 
magnitude correction method, magnitudes of 
specimen vectors were adjusted using the 
variables MXH and MXW following the cri- 
teria presented for M. audaculus, and only 
the results using this method of size correc- 
tion are presented here. The number of spec- 


imens measured for each stratigraphic hori- 
zon, and the position of these horizons in the 
stratigraphic section, is shown in figure 2. 


ANALYZING RATES OF EVOLUTION 


Individual morphological variables were 
analyzed using the technique developed by 
Haldane (1949) and treated extensively by 
Gingerich (1985), Stanley (1985), and Stanley 
and Yang (1987). This technique treats 
changes in nonmeristic morphological char- 
acters as being essentially exponential in na- 
ture. Magnitude of change in a variable over 
time is quantified using darwins. These units 
represent change in a feature by factors of e 
(~2.7) per million years. Rates are calculated 
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using the formula (In x, — In x,)/t, where x, 
and x, are initial and final mean morpholo- 
gies for a single variable and ¢ is time in mil- 
lions of years (Stanley and Yang, 1987). Both 
raw and size-corrected data were analyzed, 
in recognition that most of the phyletic change 
that occurs within species lineages is actually 
size change (Hallam, 1978; Stanley, 1985; 
Stanley and Yang, 1987). 

Because of the paucity of radiometric dates 
aside from those bracketing the upper and 
lower horizons of the Hamilton Group, anal- 
yses of rates of change in individual shell 
variables used the following temporal dura- 
tions: each formation persisted for an equal 
amount of time (1.25 Ma), and formations 
were divided up into members of equal du- 
ration. These durations are roughly support- 
ed by data presented in Brett and Baird (1986). 

Analyses of rates were also used to relate 
phenotypic evolution to intensities of selec- 
tion pressure. Although it is perhaps impos- 
sible to discern the actual potency of selective 
forces acting at the organismal level in the 
fossil record, Lande (1976) has developed a 
method for estimating the minimum selec- 
tive mortality that would produce the ob- 
served rate of phenotypic evolution. This 
method should not be viewed as a way of 
reconstructing actual selection pressures in 
the fossil record. However, it does provide a 
means of quantitatively relating patterns of 
phenotypic change through time to a hypo- 
thetical selection pressure, in order to assess 
whether those intensities of selection alone 
are what produced the given trend. If nearly 
infinitesimal selection pressures are indicat- 
ed, it suggests that less importance should be 
ascribed to such processes. 

Lande (1976) derived the following for- 
mula for estimating the proportion of the 
population culled each generation: 


EV —2 In(V 22|z|/c) 
b = ——________ (1) 
h?t 

The proportion is obtained using the estimate 
band tables of integrals for the standard nor- 
mal distribution. |z| is the difference be- 
tween In(x.) — In(x,) where x, and x, are the 
means of a population for a single variable 
at times 2 and 1. a is the mean of the standard 
deviation of each variable from those time 
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Fig. 4. Variables taken from the pedicle valve 
of specimens of Mediospirifer audaculus (Conrad). 
A, Ventral view: H, height of shell or distance from 
umbone to midline of fold; W, total width of shell; 
IAW, distance between the two ribs on each side 
of fold. B, Posterior view: HTR, orthogonal dis- 
tance between umbone and brachial valve, span- 
ning the interarea; WTR, distance from umbone 
to brachial valve, following trace of margin be- 
tween delthyrium and interarea. C, Anterior view: 
MAXF__SW, width of fold; MAXF__SH, depth 
of fold at base of shell. See text for J. 
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Fig. 5. Measurements taken from the pedicle 
valve of specimens of Athyris spiriferoides (Eaton). 
A, Ventral view: MXH, height from umbone to 
medial portion of fold; MXW, maximum width 
of pedicle valve. B, Anterior view: MFW, width 
of fold; MFH, depth of fold. See text for H and 
WU. 


intervals, h? is the heritability of the char- 
acter, and ¢ is the time in generations over 
the interval of the study. Populations must 
be checked for homogeneity of variances 
(Lande, 1976). This requirement was not sat- 
isfied for all populations. 

Values of heritability of characters (h7) are 
difficult to deduce for taxa long extinct whose 
extant brachiopod relatives have received lit- 
tle study in this area. Thus, the heritability 
values follow those of Lande (1976) and Ja- 
blonski (1987), with h? = 0.5; however, 
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changes in this variable do not strongly in- 
fluence the results because this equation 1s 
dominated by ¢, the time variable. ¢ is related 
to the number of generations over the time 
interval for which phenotypic evolution is 
being studied. Again, since establishing gen- 
eration times for these taxa is impossible, 
they are treated as having one generation per 
year. 


STATISTICAL ANALYSES 


These data were subjected to a principal 
components analysis and a canonical dis- 
criminant analysis using the PRINCOMP and 
CANDISC procedures of SAS (1987). These 
statistical analyses were conducted on sam- 
ples subdivided both by their occurrence in 
a particular stratigraphic horizon and in a 
particular paleoenvironment/biofacies. 

For the principal components analysis, only 
the variance-covariance matrix is shown. 
When that matrix is used, variables with 
higher variance are given greater weight in 
the analysis (Neff and Marcus, 1980). When 
the correlation matrix is used, the loadings 
are nearly identical, though the order of the 
first two principal components obtained is 
switched in both species. 

In order to apply the second statistical 
technique used, canonical discriminant anal- 
ysis, all groups must be multivariate normal 
and have homogeneous variance-covariance 
structure. Normality for all variables was as- 
sessed using PROC UNIVARIATE NOR- 
MAL (SAS, 1987). Although univariate nor- 
mality of all variables does not guarantee 
multivariate normality, it is a reasonable first- 
order exploratory technique for its determi- 
nation (Neff and Marcus, 1980). Across all 
horizons, variables were normally distribut- 
ed; however, within certain horizons, the dis- 
tributions of certain variables departed from 
normality. 

In addition to these two multivariate sta- 
tistical techniques, analyses using univariate 
statistics were performed. t-tests (PROC 
TTEST in SAS [1987]) were performed in 
pairwise comparisons in conjunction with the 
Bonferroni inequality to ascertain statistical- 
ly significant departures between the means 
of horizons. A value of 0.05 or less was taken 
to be statistically significant. A series of pair- 
wise comparisons were made using some of 
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the different variables for some of the key 
stratigraphic horizons. 

Additional statistical analyses were also 
conducted to take advantage of the fact that 
these brachiopods were accretionary organ- 
isms, and a history of their growth is pre- 
served in their shells. This allowed the mea- 
surement of several different ontogenetic 
stages within a particular organism. Such an 
analysis made it possible to detect any changes 
in the ontogenetic trajectory of organisms be- 
tween stratigraphic horizons and within bio- 
facies. 

Data from ontogenetic series of both bra- 
chiopod species were analyzed using an anal- 
ysis of covariance, the PROC GLM of SAS 
(1987). However, in these organisms it is im- 
possible to determine the age of an individual 
at any particular ontogenetic stage. McKin- 
ney (1986) elaborated the concept of allo- 
metric heterochrony, which uses size as a 
proxy for age/time. This allows the values of 
certain characters to be calibrated between 
individuals and across stratigraphic hori- 
zons. Ontogenetic stages which occur at the 
same size are taken to be equivalent. How- 
ever, this methodology may have problems 
since trait changes are being compared as a 
function of size instead of time, and this ig- 
nores the fact that size may be a nonlinear 
function of time (McKinney, 1988). 

In addition, Blackstone and Yund (1989) 
elegantly demonstrated that in analyses 
searching for heterochrony different results 
are often produced when either allometry or 
chronology is used to calibrate the timing of 
developmental events. Finally, Blackstone 
(1987) showed that the concept of allometry 
may not be the best way of looking at evo- 
lutionary alterations in development. There- 
fore, studies based on allometric hetero- 
chrony might not be able to illuminate 
patterns of heterochrony, though they pro- 
vide the best possible insight into these ques- 
tions, considering the type of taxa under study. 
A meaningful analysis of changes in variables 
with relative or absolute size can still be con- 
ducted, and looking at trait change as a func- 
tion of size can be of interest (Blackstone and 
Yund, 1989; McKinney, 1986, 1988). 

Because measures of absolute size are typ- 
ically used as calibrations of “‘time”’ in allo- 
metric ontogeny, this analysis uses the vari- 
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able H, the height of the shell valve, which 
correlates well with absolute size on the basis 
of the analysis presented in the section on 
MEASUREMENTS TAKEN ..., to deter- 
mine the “‘age”’ of an organism at a particular 
growth stage. Three different ontogenetic 
stages were measured for every specimen. At 
each of these three stages, four measure- 
ments, including the measurement calibrat- 
ing time, were taken. Measurements were 
taken at roughly one-quarter, one-half, and 
three-quarters of the shell height to avoid in- 
troducing possible biases. Of course, as we 
are using size to calibrate time, these data 
were not corrected for size. 

Data from ontogenetic series were sub- 
jected to an analysis of covariance. This anal- 
ysis used time as the independent variable, 
and plotted values of the dependent vari- 
ables, or morphometric measurements, 
against time. A regression line was fit to this 
scatter of data. The goal of this analysis was 
to compare the changes in variables with age 
or “time’’ across different stratigraphic ho- 
rizons and different biofacies. Differences in 
the average size of organisms across strati- 
graphic samples could lead to differences in 
the slope of regression lines between hori- 
zons, if measurements captured different por- 
tions of a nonlinear allometric growth curve. 
However, this is not a problem, because for 
all horizons the scatter of points is linear. 

For Mediospirifer audaculus, the ontoge- 
netic trajectory of three measurements was 
calculated with respect to time (H). These 
measurements are W, the width along the 
interarea from the point where the growth 
line hits the boundary between the ribs, 14 W, 
the width of the fold at the particular stage 
of ontogeny, and /, the distance from the point 
where the growth line hits the boundary be- 
tween the ribs and the interarea to the del- 
thyrium, following the trace of the growth 
line (see fig. 4A—C). 

An analysis of covariance was also per- 
formed for three variables plus the size or 
time variable (H) taken on specimens of 
Athyris spiriferoides. The variables used were: 
W, the distance between the points where the 
ontogenetic marker or growth line hits the 
margins of the pedicle valve; WU, the dis- 
tance from the position of the growth line at 
the margin to the umbone; and HTOW, the 
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TABLE 1 
Principal components analysis using the variance-covariance matrix and data corrected for size show- 
ing the eigenvalues (EV) and the proportion of the variance in the data the eigenvectors explain (PROP) 
for (A) Mediospirifer audaculus (Conrad) with only the first three eigenvectors shown, and (B) Athyris 
spiriferoides (Eaton) with only the first two eigenvectors shown. 


(A) EV PROP H IAW HTR WTR MAXF__SW MAXF__SH 
PRIN 1 .0180 .69 —.05 .O5 .69 72 05 05 
PRIN2 .0037 ad =. 59 34 — .06 —.07 68 34 
PRIN3 .0031 12 A —.27 -.01 —.01 —.17 92 
(B) EV PROP MXH MXW MFW MFH 

PRIN] .003 49 TS i, 87 .46 

PRIN2 .002 33 mers, 61 29 =.15 


distance between the point where the growth 
line was calibrated using the measurement H 
to the point where the growth line hits the 
margin of the pedicle valve (see fig. 5A,B; 
HTOW not shown). 

All data, and the SAS (1987) programs 
written to analyze these data, are available 
from the first author. 


DEFINING STASIS AND 
RECOGNIZING TAXA 


The test for stasis in species lineages is re- 
stricted to morphological changes accrued 
from earliest to latest Hamilton times (the 
stratigraphic ranges of both taxa). If there is 
no net discernible morphological change be- 
tween the end members from the bottom and 
top of the section, then the hypothesis of sta- 
sis is not rejected. Change in a species lineage 
must be irreversible to count as true evolu- 
tionary change (see Gould and Eldredge, 
1977), and oscillations are held to be equiv- 
alent to net evolutionary stability (Wnght, 
1931; Lande, 1986). This study is predicated 
on the search for anagenetic change. We do 
not claim to be able to evaluate Lande’s (1986) 
or Eldredge and Gould’s (1972) models of 
speciation. 

Because this study focuses on species-level 
taxa, it is crucial that we explicitly define a 
species. We are using the epistemological spe- 
cies concept of Eldredge and Cracraft (1980), 
which defines a species as the smallest diag- 
nosably distinct cluster of putatively repro- 
ductively interacting organisms that can be 
recognized and defined by at least one apo- 
morphy. For both species, several characters 


are available to make hypotheses asserting 
their monophyly. The spiriferid brachiopod 
Mediospirifer audaculus is diagnosed by: large 
sweeping growth lines of low convexity, large 
and extensive interarea of the pedicle valve 
with horizontal striae, large hinge teeth de- 
veloped as a steep wall, and its micro-orna- 
ment which features widely spaced growth 
lamellae. Athyris spiriferoides is diagnosed by: 
pedicle valve with small pointed teeth with 
swelling at their dorsal ends, weak dental 
plates diverging at about a 45° angle, and 
weak, oval adductor muscle scars. None of 
the characters used to define taxa were ana- 
lyzed to avoid logical inconsistency following 
Lande (1986), for if these characters changed 
they would no longer indicate the same spe- 
cies. 


RESULTS 


ANALYSIS OF TERMINAL ONTOGENETIC 
STAGES SUBDIVIDED BY 
STRATIGRAPHIC HORIZON 


Mediospirifer audaculus. — The first ex- 
ploratory technique used was principal com- 
ponents analysis. The first three principal 
components account for 95% of the variance 
in the data (see table 1A). The first principal 
component has large loadings on HTR and 
WTR, indicating that a large portion of the 
variance in the data set can be attributed to 
differences in these variables. The other vari- 
ables have near zero loadings. The second 
principal component is dominated by 
MAXF__SW and H, and the third by 
MAXF__SH. 

A canonical discriminant analysis with the 
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Fig. 6. Results from canonical discriminant 
analysis of Mediospirifer audaculus that employed 
stratigraphic horizons to discriminate specimens. 


LIEBERMAN ET AL.: STASIS IN TWO BRACHIOPOD SPECIES 11 


stratigraphic horizons used as the classes or 
groups indicates that specimens from the 
lower and uppermost portion of the Hamil- 
ton Group are very similar (see fig. 6A, table 
2A), indicating overall net stasis. However, 
changes in morphology occur within the sec- 
tion. The greatest differences are between the 
lowest stratigraphic sections, A, B, and C, and 
the intermediate sections, E, F, and G, with 
little overlap of the specimens from the Mar- 
cellus (horizon A) and Wanakah (horizon E) 
formations (see table 2A and fig. 6B). Thus, 
specimens of M. audaculus from the oldest 
sections, the Oatka Creek and Mount Marion 
shales in the Marcellus Formation, bear a 
greater resemblance to specimens from the 
youngest section, the Windom Member of 
the Moscow Formation than to those from 
the intermediate Wanakah Shale of the Lud- 
lowville Formation. 

There is considerable overlap between 
samples from successive horizons except in 
two cases. The first is between C and D, which 
are quite distinct, with a Mahalanobis dis- 
tance (D?) of 4.65. However, samples from 
the 100 m of rock between sections C and D 
could not be obtained. One might expect a 
substantial morphological gap between these 
sections due to the large interval of time for 
which no M. audaculus specimens were an- 
alyzed. The other large morphological gap 
separates samples H and J. In this case there 
is no major “sampling unconformity” sepa- 
rating the two horizons. 

The similarity between lower Hamilton and 
upper Hamilton M. audaculus, and the di- 
vergence of intermediate Hamilton M. au- 
daculus, is chiefly governed by the pattern of 
shifts in the size of the interarea. However, 
it is impossible to ascribe a cutoff point de- 
marcating those specimens with larger inter- 
areas from those with smaller interareas be- 
cause there is continuous variation between 


— 


Shown are canonical variate axes 1 and 2 and the 
scores of specimens from stratigraphic horizons: 
A, from the lowest and uppermost Hamilton 
Group, stratigraphic horizons A and K. Notice the 
great degree of overlap, and B, stratigraphic ho- 
rizons A and E. Notice the prominent separation 
between the specimens from these two horizons. 
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TABLE 2 
(A) Canonical discriminant analysis for M. audaculus showing the Mahalanobis distances (D”) between 
the different stratigraphic horizons after all vectors were corrected for size. (B) The total canonical 
structure for canonical variates 1-3 and all measured variables, with the proportion of the data explained 
by each canonical variate designated as PRO. 


(A) B Cc D E F G H J K 
A 0.34 2.24 3.11 6.17 5.77 4.54 3.77 2.43 2.19 
B 0 1.28 2.04 4.39 4.08 2.84 2.89 1.04 0.95 
Cc 0 4.65 6.74 5.65 3.40 5.64 0.83 1.33 
D 0 0.76 0.80 0.99 0.34 2.80 1.75 
E 0 0.58 1.31 0.72 4.09 2.53 
_*-F 0 0.42 0.73 3.78 2.59 
G 0) 1.23 2.26 1.67 
H 0 3.76 2.44 
J 0 0.63 
(B) PRO H IAW HTR WTR MAXF_SW MAXF__SH 
CAN1 0.71 —0.35 0.31 0.94 0.90 0.29 0.16 
CAN2 0.19 —0.08 0.59 —0.19 —0.26 —0.17 —0.01 
CAN3 0.09 0.31 —0.16 —0.01 —0.02 —0.34 0.85 


the two end points. The structure of the ca- 
nonical variates indicates that approximately 
70% of the differences of centroids for strati- 
graphic horizons can be explained by the first 
canonical variate, and HTR and WTR have 
very large loadings on this axis. The second 
canonical variate summarizes 19% of the 
variance among centroids of horizons and is 
controlled mainly by the variable [4 W (see 
table 2B). 

For the purposes of brevity and clarity, not 
all of the results obtained from comparisons 
between stratigraphic horizons using uni- 
variate statistics are presented here. How- 
ever, for those horizons found to differ by 
substantial D? values, t-tests demonstrate that 
the differences in the mean value of H7TR and 
WTR for horizons A, B, C and E, F, and G 
were Statistically significant differences in- 
dividually at the 0.0001 level, and thus using 
the Bonferroni inequality, the differences were 
significant at the 0.001 level. For those ho- 
rizons differing by small D? values, differ- 
ences in the mean value of H7R for horizons 
A and B were not statistically significant, and 
differences between horizons A and K for 
ATR were also not significant. Nonparamet- 
ric tests mirrored these results. 

Thus, statistical analysis of M4. audaculus 
indicates net stasis with some oscillation in 
morphology. The differences in HTR and 


WTR seem to be controlling much of the 
change seen in this species lineage. This pat- 
tern may be due to incipient sympatric spe- 
ciation or it may represent dynamics between 
two closely related forms. However, because 
there seems to be continuous variation in the 
size of the interarea, M. audaculus was not 
split into two separate lineages. If the pattern 
seen is within one species then we may be 
seeing changes within different populations. 
The specimens from the intermediate layers 
may reflect the increased presence of con- 
specifics from a different population that 
evolved elsewhere and migrated in, only to 
migrate elsewhere or go extinct in the Mos- 
cow Formation. As a final possibility, these 
oscillatory changes may just represent oscil- 
lation in morphology producing no net trend 
and caused by drift or wavering selection 
pressures. 

Athyris spiriferoides. — For the first principal 
component derived using the variance-covari- 
ance matrix MFW constitutes the primary 
loading (see table 1B). The second principal 
component is most affected by MXH and 
MXW. Together these axes explain 82% of the 
variance in the data. 

Canonical discriminant analysis revealed 
overall net stasis, and specimens from the 
basal and upper Hamilton Group still over- 
lap in the scores on the canonical axes (see 
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Fig. 7. Results of canonical discriminant anal- 
ysis of Athyris spiriferoides that employed strati- 
graphic horizons to discriminate specimens. Shown 
are canonical variate axes 1 and 2 and the scores 
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fig. 7A). The oldest samples, collected from 
the basal layer of the Oatka Creek Shale, dif- 
fered most from the samples in horizons O, 
P, and Q (see table 3 and fig. 7B). These sam- 
ples represent the Jaycox Shale, the Deep Run 
Member, and the Kashong Shale, respective- 
ly. For this species, unlike M. audaculus, there 
is considerable overlap between Marcellus 
formation samples and those from the Wana- 
kah Shale. Between successive samples there 
is considerable overlap, and even samples 
separated by hundreds of feet of rock, i.e., 
upper Marcellus (B) and Centerfield (C), are 
quite simular. 

Canonical variate structure suggests that 
height and negative width predominate on 
axis | and MFW and MFH predominate on 
axis 2 (see table 3B). 

t-tests on normally distributed data indi- 
cate differences in the variables MXH and 
MFW between the most divergent horizons 
Aand O significant individually at the 0.0001 
level of confidence, which, when used in con- 
junction with the Bonferroni inequality, im- 
ply differences significant at the 0.02 level. 
When the Bonferroni inequality is applied to 
differences for XW for these horizons, they 
are found to be significant at the 0.07 level, 
i.e., not significant. The nonparametric Wil- 
coxon test revealed differences between ho- 
rizons A and O for MFH only significant at 
the 0.81 level. For earliest and latest hori- 
zons, A and §, f-tests in conjunction with the 
Bonferroni inequality indicated that the mean 
values of MXH, MXW, MFW, and MFH did 
not significantly differ. 


QUANTIFYING CHANGES IN 
INDIVIDUAL VARIABLES 


Rates of change in individual morpholog- 
ical variables were calculated for both the 
total duration of the species and between the 


_ 


of specimens from stratigraphic horizons: A, A 
and S. A and S are the lowermost and uppermost 
horizons, respectively; notice their significant 
overlap. B, A and Q. Notice the prominent sep- 
aration between specimens from these two hori- 
zons. 


14 AMERICAN MUSEUM NOVITATES NO. 3114 


TABLE 3 
(A) Canonical discriminant analysis for A. spiriferoides showing the Mahalanobis distances (D?) be- 
tween the different stratigraphic horizons after all vectors were corrected for size. (B) The total canonical 
structure for canonical variates 1 and 2, with the proportion of the data explained by each canonical 
variate designated as PRO. 


(A) B Cc D E F G H I 
A 0.58 3.84 4.66 0.25 2.86 4.47 1.97 1.61 
B 0 2.08 2.26 0.41 1.58 2.88 0.95 1.26 
C 0 0.28 3.16 0.41 0.78 0.73 1.14 
D 0 3.76 0.86 1.21 1.18 2.09 
E 0) 1.94 3.13 1.23 1.09 
F ) 0.23 0.15 0.38 
G 0 0.63 0.93 
H 0 0.24 

K L M N oO P Q R S 
A 207 3.90 3.54 4.03 6.22 5.76 7.88 2.47 3.09 
B 1.96 1.98 2.23 3.31 3.75 3.99 5.50 1.05 2.07 
C 1.44 0.24 0.56 1.39 0.45 1.48 1.43 0.43 1.07 
D 2.48 0.36 1.05 2.43 0.53 1.98 1.82 0.67 1.75 
E 1.51 3.39 2.48 2.96 5.32 5.84 7.53 1.90 2.13 
F 0.50 0.80 0.20 0.63 1.17 2.69 2.83 0.37 0.27 
G 0.77 1.40 0.24 0.49 1.27 3.60 3.23 1.02 0.43 
H 0.59 0.88 0.61 1.10 1.60 2.69 3.25 0.22 0.24 
I 0.19 1.53 0.73 0.74 2.38 3.16 3.87 0.72 0.39 
K ) 2.24 0.50 0.27 2.89 4.33 4.74 1.35 0.62 
L 0 1.36 2.46 0.34 0.78 0.99 0.27 1.27 
M ) 0.34 1.57 3.62 3.46 1.01 0.78 
N 0 2.55 4.62 4.41 1.87 0.87 
O 0 1.07 0.61 0.93 1.58 
P 0 0.36 1.54 2.91 
Q 0 2.08 3.81 
R 0 0.60 
(B) EV PRO MXH MXW MFW MFH 
CANI 0.48 0.62 0.96 —0.95 —0.20 0.17 
CAN2 0.21 0.26 0.22 -0.21 0.98 0.45 


most divergent horizons, determined using 
Mahalanobis D? values, to recover the max- 
imum possible rates of change. 
Mediospirifer audaculus. — The results 
from raw data, which include evolutionary 
changes in size, are presented in table 4A, 
and indicate modest levels of change. A 
change of 2.7% per million years represents 
10 millidarwins (Stanley and Yang, 1987). 
These values are comparable with those for 
changes in size discovered by Hallam (1975) 
from his studies on Jurassic bivalve lineages. 
However, when the data are corrected for size 
the amount of change measured in millidar- 
wins falls precipitously (see table 4A) such 


that overall changes are very small, and com- 
parable to the minimal changes recovered by 
Stanley and Yang (1987) from their analysis 
of changes in Neogene bivalves. Thus, much 
of the change in these data is related to changes 
in size. Still, it is important to point out that 
when only earliest and intermediate sections 
are considered, significant changes do occur 
in the variables HTR and WTR, but these 
changes are reversible. 

Athyris spiriferoides. — The results from 
raw and size-corrected data are presented in 
table 4B. Again, very little change transpires 
from earliest to latest occurrence of the spe- 
cies. 
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TABLE 4 

Rates of change in individual morphological 
variables in millidarwins (MD) between earliest 
and latest horizons in the Hamilton Group and 
also between those horizons which differed the 
most on the basis of the Canonical Discriminant 
Analysis. (A) Mediospirifer audaculus, raw data 
and data corrected for differences in size. (B) Athyris 
spiriferoides, raw data and data corrected for dif- 
ferences in size. 


MD of change between horizons 


(A) Raw data Corrected data 

Variables A-K A-E A-K A-E 
H 85 167 2 9 
[IAW 103 210 14 32 
HTR 100 309 17 127 
WTR 92 290 9 111 
MAXF__SW 87 184 4 6 
MAXF__SH 120 240 30 65 
(B) MD of change between horizons 

: Raw data Corrected data 

Vari- ee ee ee a a 

ables A-S A-O A-Q AS A-O A-Q 
MXH 87 100 134 2 21 23 
MXW 89 126 170 ns 15 15 


MFW 143 = 209 230 55 94 78 
MFH 79 144 165 15 29 5 


ASSESSING RELATIONSHIP BETWEEN 
ENVIRONMENT AND STASIS 


To assess the role environment may play 
in mediating stasis or change in species mor- 
phology, we concentrated on morphological 
changes through time within a single paleoen- 
vironment/biofacies. This was done to as- 
certain if there are any morphological shifts 
through time in any biofacies. If constant en- 
vironment exerts constant selection pres- 
sures, and thus helps mediate stasis via sta- 
bilizing selection, then one would predict that 
through time there would be little morpho- 
logical change in a species, as long as the 
environment remained constant. Canonical 
discriminant analyses were used to assess 
morphological differences between strati- 
graphic horizons within a particular biofa- 
cies. 

Mediospirifer audaculus. — Analyses of 
morphological change within single paleoen- 
vironments/biofacies were conducted for 
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TABLE 5 

Canonical discriminant analysis for M4. auda- 
culus showing the Mahalanobis distances (D?) (A) 
between different stratigraphic horizons for pa- 
leoenvironment/biofacies B, (B) the distances be- 
tween different stratigraphic horizons for paleoen- 
vironment/biofacies D, and (C) between the 
different paleoenvironments/biofacies after all 
vectors were corrected for size. 


(A) D E F H 
Cc 9.76 13.28 8.91 7.26 
D 0 1.45 5.54 3.13 
E 0 FQ. 6.02 
F 0 5.29 
(B) F G H J K 
E 0.62 4.13 2.71 14.59 8.64 
F 0 2.52 1.66 12.24 6.03 
G 0 2.09 5.89 1.32 
H 0 9.93 2.98 
J 0 5.59 
(C) Cc D E F 
B 1.80 0.56 0.58 1.61 
Cc 0 3.21 1.23 2.21 
D 0 0.82 2.02 
E 0 0.91 


biofacies B, C, and D. (The other biofacies 
were present in less than three of the strati- 
graphic horizons considered herein, and thus 
could not be subjected to a canonical dis- 
criminant analysis.) There are significant 
shifts in morphology within all of these bio- 
facies through time (see table 5A, B). For 
example, consider the Mahalanobis distance 
between stratigraphic horizons E and J from 
all biofacies (table 2A, fig. 8A), and the dis- 
tance between stratigraphic horizons E and J 
solely for specimens from biofacies D (table 
5B, fig. 8B). There has been more change 
through time within a single paleoenviron- 
ment/biofacies than across all biofacies. These 
results are reinforced when patterns of change 
in individual variables are examined. For in- 
stance, the changes in all variables, particu- 
larly HTR and WTR, that transpired between 
the interval of time in which horizon E and 
horizon J was deposited (approximately 1.5 
Ma) is smaller when all paleoenvironments/ 
biofacies are considered than when only a 
single biofacies is considered (see table 6A, 
B). 
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Athyris spiriferoides. — Analyses to sample 
temporal variation again show that there are 
larger differences manifested within a single 


biofacies through time than when all biofa- 
cies are considered together (see table 7A). 
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Fig. 8. Results of canonical discriminant anal- 
ysis of specimens of Mediospirifer audaculus. 
Shown are canonical variate axes | and 2. A, B, 
Using stratigraphic horizons to discriminate spec- 
imens. A, Specimens from all environments were 
considered. The scores of samples from horizons 
E and J are shown. B, Only specimens from bio- 
facies D were considered, scores of specimens from 
horizons E and J shown, compare with fig. 8A; 
notice the greater difference between the two ho- 
rizons here. C, Using occurrence in a particular 
paleoenvironment/biofacies to discriminate spec- 
imens. Shown are the scores of specimens from 
paleoenvironments B and D on these axes. Notice 
the near complete overlap of specimens from these 
two paleoenvironments. 


Analyses of changes in individual variables 
are comparable to those from M. audaculus 
and are not presented here for the purposes 
of brevity and clarity. 

These results from both taxa suggest that 
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TABLE 6 

Rates of change in individual morphological 
variables in millidarwins (MD) for Mediospirifer 
audaculus between stratigraphic horizons E, the 
lower Wanakah Shale, and J, the lower Windom 
Shale, a 1.5 Ma interval, (A) when all paleoen- 
vironments/biofacies are considered and (B) when 
only patterns of change in the Pseudoatrypa pa- 
leoenvironment/biofacies (biofacies D) are con- 
sidered. 


MD of MD of 
(A) change (B) change 
H 0 H 7 
IAW 30 TAW 0 
HTR 46 HTR 406 
WTR 40 WTR 560 
MAXF__SW 7 MAXF_SW 20 
MAXF__SH 53 MAXF_SH 133 


the environment alone may not play a prom- 
inent role in maintaining morphological sta- 
bility. In each case we are seeing moderate 
morphological changes through time within 
paleoenvironments/biofacies. However, 
when all biofacies are lumped together, we 
see little directed net change. This argues 
against the notion that morphologically sta- 
ble populations occupy the same environ- 
ment through time, and therefore questions 
the exclusive role that stabilizing selection 
plays in mediating stasis. It suggests that, for 
the limits of resolution of this analysis, in 
these taxa the argument for stasis elaborated 
by Stanley (1979) and Eldredge (1989) and 
based on Wright’s (1931) view of species is 
supported. Wright (1931) recognized that 
species are generally broken up into different 
groups or demes that live in different envi- 
ronments. Stanley (1979) and Eldredge (1989) 
suggested that each of these demes will un- 
dergo largely independent adaptive histories 
in reaction to local environmental conditions 
so that the net sum of these different histories 
will generally lead to no net change. 

We do not claim to be analyzing actually 
reproducing populations or demes in the 
Hamilton Group. However, we are consid- 
ering statistical populations of organisms oc- 
curring in similar environments. We refer to 
these samples as environmental populations 
(see DISCUSSION). In the Hamilton Group 
these environmental populations appear to 
undergo significant changes in morphology 
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TABLE 7 
Canonical discriminant analysis for A. spirifer- 
oides showing the Mahalanobis distances (D”) after 
all specimens were corrected for size between sam- 
ples from (A) the different stratigraphic horizons 
from paleoenvironment/biofaces B and (B) differ- 
ent paleoenvironments/biofacies. 


(A) E F H K 

C 4.42 2.12 3.58 3.69 

E 0 2.03 2.20 1.47 

F 0 0.77 0.44 

H 0 1.16 

(B) Cc D E F G 
B 0.22 0.25 2.97 0.55 1.25 
Cc 0 0.48 1.95 0.64 0.65 
D 0 2.43 0.12 0.85 
E 0 2.07 0.65 
F 0 0.57 


through time. However, the net sum of these 
morphological changes across all environ- 
mental populations leads to no net shift in 
species morphology. 

No prominent geological features have been 
recognized in the Middle Devonian of New 
York State that would have divided the rang- 
es of the species being studied into different 
geographic subregions. Therefore, this vari- 
ation may be independent of geographic iso- 
lation. However, different environments were 
in relatively different positions on the sloping 
profile of the shoreline shown in figure 3. 

These differences in patterns of phenotypic 
evolution in single paleoenvironments and 
across all paleoenvironments can be quan- 
tified using Lande’s (1976) procedure for cal- 
culating minimum selective mortalities using 
equation (1) (see table 8A, B). The values of 
b correspond to the truncation point for mor- 
phologies, scaled as the number of standard 
deviations from the average phenotype. The 
number of selective deaths per individual per 
generation is also given. This is shown only 
for M. audaculus here, but a similar picture 
emerges for A. spiriferoides. Minimum selec- 
tive mortalities were calculated for the mor- 
phological variables that differed most be- 
tween horizons E and J. Mortalities are 
presented both across all paleoenvironment/ 
biofacies and only in biofacies D. These re- 
sults indicate that selective mortality in re- 
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TABLE 8 

Differences in patterns of phenotypic evolution 
in Mediospirifer audaculus between stratigraphic 
horizons E, the lower Wanakah Shale, and J, the 
lower Windom Shale related to minimum selec- 
tive mortalities using equation (1) in the text, (A) 
across all paleoenvironments/biofacies, and (B) 
only for the Pseudoatrypa biofacies (biofacies D). 
Shown are values for |z|, ¢, b, and approximate 
minimum selective mortalities (MSM), all defined 
in the text. Only the most divergent morphological 
variables are presented. 


(A) \z| o b MSM rates 
HTR 0.07 0.06 4.99 3x 10-7 
WTR 0.06 0.06 5.02 3x 10-7 
MAXF_SH 0.08 0.06 4.96 3 x 10-7 
(B) \z| o b MSM rates 
HTR 0.61 0.046 4.48 3.4 x 10-6 
WTR 0.56 0.047 4.50 3.4 x 10-6 
MAXF__SH 0.08 0.065 4.80 8 x 10-7 


lation to the characters studied is extremely 
low in this species lineage, and this level of 
selection could probably not be distinguished 
from drift. For instance, across all paleoen- 
vironments considered, only three selective 
deaths per 10 million organisms occurred. 
These values are comparable to those recov- 
ered by Lande (1976) from a series of other 
fossil taxa, and suggest that across all pa- 
leoenvironments only those phenotypes 5 SD 
from the mean would be eliminated. As Hal- 
dane (1949) recognized, fossil specimens that 
differed this much from mean morphology 
would probably be assigned to a separate spe- 
cies. 

It is interesting to note that selective mor- 
tality is a factor of 10 higher in a single pa- 
leoenvironment/biofacies than across all pa- 
leoenvironments/biofacies. These selective 
mortalities are still very low, but this infor- 
mation argues against the conclusion that sta- 
bilizing selection is preserving stasis. If any- 
thing, selection pressures are stronger within 
a single environment than across all envi- 
ronments. 

It is important to note that in these anal- 
yses, environments were constant only for the 
range of variables we recognized. There are 
several other variables that we could not or 
did not measure but which may have varied. 
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ASSESSING ECOPHENOTYPIC VARIATION 


In order to assess how ecophenotypic ef- 
fects may have influenced morphology, a ca- 
nonical discriminant analysis was used to 
partition specimens according to the paleoen- 
vironment/biofacies they occurred in. If eco- 
phenotypic effects strongly influenced species 
morphology in these brachiopod taxa, then 
the scores of specimens on canonical variates 
axes should significantly differ if specimens 
hail from different biofacies. However, for 
both M. audaculus and A. spiriferoides the 
opposite result was found. Scores of speci- 
mens from different paleoenvironments/bio- 
facies completely overlapped on canonical 
variate axes (fig. 8C), and there was little dif- 
ference between biofacies (tables 5C, 7B). 
Thus ecophenotypic effects do not play a 
prominent role in these data. 


STATISTICAL ANALYSIS OF 
ONTOGENETIC TRAJECTORIES © 


Mediospirifer audaculus. — An analysis of 
covariance, when conducted to test for het- 
erogeneities of slope in plots of the variable 
W versus time (#) across different strati- 
graphic horizons gave the slopes in table 9A. 
When regression lines are fit to specimens 
from different stratigraphic horizons, slopes 
from the uppermost and lowermost sections 
do not differ significantly, indicating stasis in 
allometric heterochrony. However, sections 
E and H do differ significantly from both ho- 
rizons A and B and J and K (see fig. 9A). 
These differences in slope, although statisti- 
cally significant, are not large. For example, 
according to table 9A, the slope fitted to the 
regression line for stratigraphic horizon K is 
1.79 + SE. For horizon H it is 2.16 + SE. 
(No units are given for these slopes because 
both axes depict distance measurements.) 

Some statistically significant differences are 
also found for the y-intercepts of these dif- 
ferent regression lines. Stratigraphic horizons 
that show significant differences in intercepts 
from the calibrating, uppermost horizon, are 
not always the same as those which tend to 
differ in slope. 

When the covariance analysis was con- 
ducted for the variables [4 W and J, patterns 
similar to those discovered for the variable 
W emerged, but the heterogeneities of slope 
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were of lesser magnitude. For the purposes 
of brevity and clarity these results are not 
presented here. 

The results of the analysis of covariance 
for ontogenies of M. audaculus suggest two 
broad patterns: 1) there is concordance be- 
tween the allometric ontogenies of specimens 
from the lowermost and uppermost strati- 
graphic horizons, with varying differences in 
samples from the middle horizons; 2) a 15% 
difference in slopes is found between those 
samples that significantly differ. Thus, some 
oscillatory change does occur in the allome- 
tric ontogeny of this species. However, no net 
change is accrued over earliest to latest Ham- 
ilton time (5 million years). The concordance 
between uppermost and lowermost strati- 
graphic samples with differences emerging in 
the middle sections roughly mirrors the re- 
sults from the canonical discriminant anal- 
ysis. 

Athyris spiriferoides. — Analysis of co- 
variance to test for heterogeneities of slope 
in plots of the variable W versus time (#7) 
across different stratigraphic horizons indi- 
cates significant differences in slope between 
the lowermost and uppermost horizons, A 
and S. In addition, horizons E, P, and Q differ 
significantly from S. Horizons P and Q differ 
most significantly from horizon A (see table 
9B and fig. 9B). There is a significant shift in 
allometric ontogeny through time in A. spi- 
riferoides, with differences in slope of about 
20% between uppermost and lowermost sec- 
tions. However, oscillation in this trend ap- 
pears to be present, with a subtle reversion 
back to the earlier pattern of allometric het- 
erochrony in the latest horizons after 30% 
differences in slope were manifested between 
horizon A and horizons P and Q. 

The patterns retrieved when ontogenetic 
trajectories for the variables WU and HW 
were considered do not differ from those 
found for the variable W, except that the 
magnitudes of the greatest differences in slope 
are smaller, only 15%. For the purposes of 
brevity and clarity these results are not pre- 
sented here. 


DISCUSSION 


Over 1000 specimens of two species of bra- 
chiopods from Middle Devonian rocks of 
New York State were subjected to morpho- 
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TABLE 9 

The slopes obtained by an analysis of covariance 
used to assess allometric heterochrony in the 
Hamilton Group. This analysis was performed to 
test for heterogeneity of slopes for values of the 
variable W vs. time (H) from different stratigraph- 
ic horizons for (A) specimens of M. audaculus, 
and (B) for specimens of A. spiriferoides. The sig- 
nificance of the departures in slope relative to the 
reference sections (stratigraphic horizons K and S 
respectively) is shown along with the standard er- 
ror of estimate of the slope. 


(A) Slope Signif SE 
A 1.81 0.88 0.15 
B 1.72 0.72 0.21 
C 2.03 0.14 0.16 
D 1.87 0.47 0.11 
E 2.05 0.01 0.10 
F 1.96 0.07 0.09 
G 1.99 0.07 0.11 
H 2.16 0.003 0.12 
J 1.84 0.66 0.10 
K 1.79 0.07 
(B) Slope Signif SE 
A 0.99 0.001 0.08 
B 1.13 0.11 0.07 
Cc 1.23 0.70 0.06 
D 1.27 0.75 0.08 
E 1.13 0.04 0.06 
F 1.23 0.72 0.06 
G 1.22 0.51 0.06 
H 1.21 0.42 0.06 
I 1.19 0.25 0.06 
K 1.16 0.08 0.05 
L 1.27 0.81 0.06 
M 1.24 0.74 0.06 
N 1.18 0.20 0.06 
O 1.32 0.27 0.06 
P 1.39 0.01 0.06 
Q 1.40 0.03 0.07 
R 1.34 0.17 0.06 
S 1.26 0.05 


metric analysis over a variety of environ- 
ments and stratigraphic horizons. Measure- 
ments were taken to capture information both 
on overall shell shape and organismal ontog- 
eny. The results from multivariate statistical 
analysis of morphometric data and rates of 
change of individual variables from Medio- 
spirifer audaculus indicate net stasis during 
its 5 million year history. However, shifts in 
morphology and allometric heterochrony did 
occur during this time. The greatest morpho- 
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Fig. 9. Scatter plot summarizing the analysis 
of covariance which tested for heterogeneity of 
slopes of the variable W versus time (1) calibrated 
allometrically, for A, Mediospirifer audaculus and 
B, Athyris spiriferoides from different stratigraphic 
horizons. A, Horizons H and K. Notice the steeper 
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logical excursions were found between sam- 
ples from the lowest and intermediate sec- 
tions of the Hamilton Group. 

The same largely holds true for Athyris spi- 
riferoides, with overall net stasis in mor- 
phology, though not in allometric hetero- 
chrony. The changes in allometric 
heterochrony in this species fits Hallam’s 
(1978) suggestion that they may display phy- 
letic trends. Oscillations do occur within the 
section, with the greatest departures in mor- 
phology and allometric heterochrony be- 
tween oldest and intermediate horizons, 
though these intermediate horizons are youn- 
ger than those that encompassed the greatest 
departures in M. audaculus. 

Analyses also indicate that ecophenotypic 
effects do not have an important influence on 
the morphology of either species. In addition, 
the morphology of both Mediospirifer au- 
daculus and Athyris spiriferoides, when quan- 
tified using multivariate morphometrics and 
rates of change of individual characters, was 
found to change more through time within a 
single paleoenvironment/biofacies than when 
all biofacies are considered, and this may have 
implications for the way we view stasis. It 
suggests that populations are not morpho- 
logically stable in the same paleoenviron- 
ment/biofacies over long intervals of time. 
Typically the environment, via stabilizing se- 
lection, has been viewed as an important force 
contributing to stasis of species lineages. 
However, these changes in the face of a con- 
stant environment (though not necessarily a 
constant selective regime, but assuming the 
two are related) suggest that environmental 
stability does not necessarily preserve species 
morphology. Amounts of phenotypic selec- 
tion on the variables studied were found to 
be a factor of ten stronger in a single paleoen- 
vironment than when all paleoenvironments 
were considered. This suggests that stabiliz- 
ing selection is not the sole factor contrib- 
uting to stasis in these taxa. 

An additional mechanism for stasis emerg- 
es from the way in which species are orga- 


, 


slope of H. B, Horizons A and P. Notice the steep- 
er slope of P. 
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nized. Wright (1931) recognized that species 
are typically broken up into reproductive 
groups occurring in separate environments, 
and these groups, referred to as demes, will 
undergo partly separate adaptive histories. 
The morphology of organisms in these demes 
may change through time as they adapt to 
local conditions. However, the net sum of 
these changes will often cancel out such that 
the overall morphology of the species re- 
mains constant (Stanley, 1979; Eldredge, 
1989). This is intuitively obvious, since only 
if all morphological changes in all environ- 
ments were in the same direction would net 
change in species morphology be significant. 

We do not claim to be studying actual 
demes or reproductive populations, in the 
restrictive sense of Wright (1931). Instead, 
we are studying statistical populations of or- 
ganisms occurring in similar paleoenviron- 
ments through time. We refer to these groups 
as “environmental populations,” and they 
represent lineages in constant environments 
in the sense of Eldredge (1989) and in the 
broad sense of Wright (1931). Study of such 
groups provides a way of assessing the rela- 
tionship between organismal morphology and 
environment. Statistical analyses suggest that 
these “environmental populations” do change 
significantly through time. These morpho- 
logical changes may be due to adaptation, 
random effects, or the invasion of popula- 
tions from unsampled geographic regions. 
However, the latter seems least likely as there 
are no prominent morphological breaks be- 
tween successive stratigraphic horizons for 
either species except in one instance for M. 
audaculus, and in this case the two horizons 
were separated by a body of rock that could 
not be sampled. 

Because a species is a unit in nature broken 
up into different subparts, the net change in 
a species is equal to the sum of the morpho- 
logical changes in all the different ‘“‘environ- 
mental populations.” Each of these parts un- 
dergoes an independent history, and thus as 
long as a species occurs in several different 
environments one would predict overall that 
it would resist change. 


CONCLUSIONS 


For the two brachiopod species studied, 
Mediospirifer audaculus and Athyris spirifer- 
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oides, the hypothesis of morphological stasis 
could not be refuted. Samples from the first 
and last occurrences of the two species, from 
the lowermost and uppermost stratigraphic 
horizons of the Hamilton Group, still show 
considerable morphological overlap. Changes 
in morphology do occur within the Hamil- 
ton, but these changes are reversible and 
therefore do not count as evolutionarily sig- 
nificant in the sense of Gould and Eldredge 
(1977) and Lande (1986). 

When single paleoenvironments/biofacies 
are studied using terminal ontogenetic stages, 
morphological changes are found to occur 
through time. These changes are larger than 
the changes that occur through time when the 
entire range of environments from which a 
species is known are considered. This implies 
that persistence in environments through time 
in conjunction with stabilizing selection is 
not the only cause of morphological stability. 
Instead, stasis may also be the result of the 
organization of species. Separate environ- 
mental populations of a species will change 
through time, but the changes in all these 
different environmental populations will gen- 
erally produce no net morphological change. 
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